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Nucleon-antinucleon annihilation in the large N c limit of QCD in the Witten regime of fixed 
velocity is considered with a focus on the spin and isospin dependence of the annihilation cross- 
section. In general, time-reversal and isospin invariance restricts the annihilation cross-section to 
depend on 6 independent energy-dependent terms. At large N c , a spin-flavor symmetry emerges in 
the theory that acts to further restrict the dependence of the annihilation cross-section to three of 
these terms; the other terms amount to 1/N C corrections. Assuming dominance of the leading order 
terms, several identities are derived that relate annihilation in different spin-isospin channels. A key 
prediction is that for unpolarized nucleons in Witten kinematics, the proton-antiproton annihilation 
cross-section should be equal to the proton-antineutron annihilation cross-section up to corrections 
of relative order 1/N C . Unpolarized nucleon-antinucleon annihilation data appears to be consistent 
with this expectation. 


I. INTRODUCTION 

Quantum chromodynamics (QCD) has long been ac¬ 
cepted as the underlying theory for hadronic physics. 
However explicit calculations for many phenomena in the 
theory remain elusive. On the one hand, perturbation 
theory in the coupling fails at small momenta. On the 
other hand, Monte Carlo evaluation of the Feynman path 
integral of a Euclidean space lattice version of QCD is 
currently the only viable non-perturbative method cur¬ 
rently available to obtain observables directly from the 
QCD Lagranagian and it has significant limitations. Lat¬ 
tice techniques are very well suited for static properties 
of hadrons. It is also possible to study low energy scat¬ 
tering of hadrons using Liischer’s methodf!®] of relating 
S-matrix information such as phase shifts to energy levels 
in a finite spatial box. While there has been considerable 
recent progress in developing this method into a practi¬ 
cal tool, it is clear that this approach becomes increasing 
unwieldy as energies increase. It is clear that the Liischer 
approach is not viable for scattering observables associ¬ 
ated with higher energy scattering for which phase space 
allows a very large number of particles in the final state. 

While a direct calculation of high-energy scattering ob¬ 
servables directly from QCD is well beyond the present 
state of the art, it is possible to learn something about 
some high-energy scattering observables from the large 
N c limit of QCD. Recall that ’t Hooft[5] recognized long 
ago that 1/IVc, where N c is the number of colors, can 
be used as an expansion parameter in QCD. The under- 
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lying premise of this approach is that a world with an 
infinite number of colors is qualitatively similar to the 
world of N c = 3. Given this premise one should replace 
the SU(3) gauge group with SU(N C ). This makes it pos¬ 
sible to study observables of interest in the limit where 
N c goes to infinity and to calculate systematic correc¬ 
tions in a power series in l/N c . The difficulty with such 
an approach is that the large N c limit of QCD, while pro¬ 
viding great simplifications compared to QCD at N c = 3, 
in general is still not tractable by analytic means except 
in special cases where gluons do not play a role as dy¬ 
namical degrees of freedom, such as the ‘t Hooft model 
in which space-time is 1+1 dimensional [6j or QCD in 
the regime in which quark masses are much larger than 
the QCD scale[7JEi- Fortunately, even when one cannot 
directly compute even the leading order behavior, it is 
often possible to deduce some properties. For example, 
it is often possible to deduce the N c scaling behavior of 
various observables which may allow for qualitative pre¬ 
dictions of the relative size of various observables—under 
the assumption that three can be considered as large for 
the purposes of the analysis. 

It turns out that among the phenomena for which pre¬ 
dictions at large N c exist are certain properties of scatter¬ 
ing observables in nucleon-nucleon scatteringj^HUS- All 
of these depend on the system being in the regime of Wit¬ 
ten kinematics [7j in which the incident momentum, along 
with the nucleon mass, scales as N )—ensuring that the 
velocity is held fixed as N c diverges. Note that Witten 
kinematics is intrinsically a high energy regime in that as 
N c goes to infinity, the initial kinetic energy grows with 
N c and the number of mesons kinematically allowed to be 
created in the scattering process also grows with N c . A 
key result is that spin-flavor dependence of a certain class 
of scattering observables jS] is constrained by an emer- 
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gent symmetry at large N c |12H19| . While the large N c 
analysis of the spin-flavor dependence of nucleon-nucleon 
scattering is more than a decade old [5], it has basically 
not been empirically tested against real world data in or¬ 
der to determine the extent to which the large N c limits 
are a reasonable caricature of the physical N c = 3 world 
with regard to scattering. There are a couple of reasons 
for this: first, to get the full spin-flavor dependencies 
one needs to have access to data for both proton-proton 
and proton-neutron scattering with polarized beams and 
targets. Such data are relatively rare. Moreover, experi¬ 
mentalists in the field do not typically analyze their data 
in terms of the types of observable for which the large N c 
analysis applies. 

There has been one attempt to analyze real world 
data[2(]| to see if the predicted spin-flavor dependence 
due to large N c can be seen at least approximately. That 
attempt exploited the fact that formally, the large N c 
analysis should apply when N c is large enough to justify 
a semi-classical treatment, which occurs when the inci¬ 
dent momenta are much larger than the typical scales of 
QCD. It is easy to see that at a formal level this condition 

can be met for sufficiently large N c , even below the elastic 
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threshold which occurs for momenta p ~ \JM n m. K ~ N£ . 
This implies that if N c were large enough, the predicted 
spin-flavor dependence should be evident in elastic scat¬ 
tering just below the elastic threshold. Empirically, the 
predicted pattern is not seen in the elastic scattering 
data even approximately [20j. This is hardly surprising. 
Firstly, in the real world in which N c = 3, it may be rea¬ 
sonable in some circumstances to argue that N c is large, 
but it stretches credulity to argue that the square root of 
three should be considered large as well. Secondly, the 
momentum scale at the elastic threshold is anomalously 

small due to the approximate chiral symmetry of QCD 
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which leads to a light pion. Thus, the value of N c 
needed to make the scattering semi-classical needs to be 
particularly large to compensate for the small size of m T . 

In this paper we extend this large N c analysis to the 
problem of nucleon-antinucleon scattering at large N c 
with an emphasis on the annihilation cross-section. In 
particular, we determine the leading order spin and flavor 
dependence of the annihilation cross-section. It is inter¬ 
esting to show that the kind of analysis used in nucleon- 
nucleon scattering is generalizable to a new class of prob¬ 
lems. More significantly, the annihilation cross-section is 
something that experimentalists conventionally measure 
and so it provides a possible testing ground for the ideas 
underlying scattering in large N c QCD. 

II. BARYONS AT LARGE N c 

The present analysis of nucleon-antinucleon annihila¬ 
tion is largely analogous to previous work on nucleon- 
nucleon scattering. It is therefore useful to review 
baryon-baryon scattering at large N c and the analysis 
of single baryons upon which it depends before turn¬ 


ing to the problem at hand. This section briefly reviews 
baryon properties at large N c and the following one re¬ 
views baryon-baryon scattering. The basic framework for 
treating both of these problems was largely laid out in 
Witten’s classic 1979 paper[7 . In that paper, intuition 
about baryons in large N c QCD is developed in the limit 
where all quarks are heavy: toq 3> Aqcd- 

A critical first piece of this analysis is the argument 
that generic baryon properties can be accurately de¬ 
scribed by a mean-field treatment. It was then argued 
that N c scaling rules such as the fact that the baryon 
mass scales with N c hold even away from the heavy quark 
limit. A key observation of Witten was that properties 
of baryons scale with N c in precisely the same way that 
properties of solitons scale with 1 /g 2 where g is the cou¬ 
pling constant. This leads to the natural expectation that 
generic properties with scaling in soliton models such as 
the Skyrme model [2TM23] faithfully reproduce the large 
N c scaling behavior of QCD, with a semi-classical treat¬ 
ment of the weakly coupled theory playing the role of 
the mean-field treatment. Following Witten’s paper it 
was rapidly realized that there were a number of pre¬ 
dictions of Skyrme type models which appeared to be 
model independent in the sense that they held regardless 
of the parameters of the models or the details of the La- 
grangian; they only depended on the symmetries of the 
theory and semi-classical treatment |24fl?7] . 

A large number of the model-independent predictions 
of baryons at large N c concern symmetries. For simplic¬ 
ity, we will restrict our attention in this paper to systems 
with two degenerate light flavors. At the level of large 
N c QCD for baryons, consistency relations imply the ex¬ 
istence of an emergent contracted SU(4) symmetry (a 
special case of the more general property of Nf degener¬ 
ate flavors yielding a contracted SU(2iV/) symmetry) jT21- 
fl9l . The algebra associated with this contracted SU(4) 
group is specified in terms of 15 generators, Ji, I a and 
X ia with i and a running from 1 to 3. The commutation 
relations are 

[Jj, Jj] = 

[IaJb\ = 

[T i , A^ a ] — 

[Ia,Xib] = 

[Jija] = 

[AC/, A jb\ 

The relevant representation of this group is infinite di¬ 
mensional and consists of states with I = J. 

The key thing which enables one to make predictions 
is that for the nucleon states of interest, the matrix el¬ 
ements of all of the generators are of order unity, that 
is Ad?. In the case of matrix elements of J l or I a this 
is obvious since the nucleons have spin and isospin 1/2. 
For the case of X ia it is less obvious but still true. It 
ultimately follows from the fact that this is a contracted 
symmetry. In effect X ia = G ia /N c where G ia ~ N c is the 
naturally arising object if there were an ordinary SU( 4) 
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rather than a contracted one. Thus, for example, the ma¬ 
trix elements of G za would be of order N c because these 
elements are (./V|'07 z 75t“'0| N). This is the strength of 
the axial vector coupling, which is of order N c due to 
normal large N c counting rules. 

Model-independent relations can then emerge when 
observables are related by Clebsch-Gordan coefficients of 
the emergent symmetry. In the Skyrme model the hedge¬ 
hog structure of a single isolated classical soliton solution 
breaks both rotational and isospin symmetry. Physically 
these hedgehogs correspond to superpositions of physical 
states with well-defined spin and isospin quantum num¬ 
bers. Since the results of interest are model-independent, 
it is sufficient to use the simplest version of the Skyrme 
model, which uses only pion degrees of freedom. These 
are collected in a single matrix-valued field U (x, t) with 
U G SU(2). The classical zero modes associated with this 
breaking become collective coordinates and must be re¬ 
quantized in order to restore the symmetries of the phys¬ 
ical states [25]. These requantized states are connected 
by the same group structure as the one in the underly¬ 
ing large N c theory and hence respect all of the same 
symmetry relations. 

A single hedgehog for a soliton with baryon number 
unity may be written as 

U(x, t) = A(t) Uo(x) T _1 (f) with (2) 

Uq(x) = e iF ^ T ' & 

where F is a radial function that satisfies F(0) = 7r, F —> 
0 as r —> oo and minimizes the energy subject to that con¬ 
straint. A{t) G SU(2) is a spatially uniform matrix which 
specifies the particular orientation of the hedgehog-that 
is, the degree to which spin and isospin are aligned. A 
can be written as A = ao + fa • r with Oq + a 2 = 1. Since 
A(t) is specified by four variables and one constraint, it 
can always be specified by three parameters- for exam¬ 
ple the Euler angles. These parameters become the col¬ 
lective coordinates for the single baryon problem. Each 
of the requantized physical states | m s rrii ) can be inter¬ 
preted as corresponding to collective wave functions in 
terms of the A. Up to overall normalization constants, 
these collective wave functions turn out to be the Wigner 
D-matrices [52]. Physical quantities are obtained from 
appropriate integrals over A weighted appropriately by 
the quantity of interest and the collective wave-functions 

Let us return to those relations that hold indepen¬ 
dently of the parameters in the model or the number 
of terms in the Lagrangian. Such relations apparently 
depend entirely on the fact that the underlying baryon 
has a hedgehog structure and the physical states are ob¬ 
tained using collective wave functions that are given by 
Wigner D- matrices [24H57] . As such, it is highly plausible 
that they are consequences of the emergent contracted 
SU(2iV/) at large N c . In numerous cases it has been 
shown that it is indeed the case that such a relation as 
is seen in the Skyrme model can be derived directly from 
the group-theoretic structure with no additional input. 


In addition there also exist many relations in Skyrme 
type models that involve the behavior of a quantity as 
the chiral limit of = 0 is approached but are other¬ 
wise independent of the parameters of the model. In all 
known cases, such behavior can be inferred from large 
N c chiral perturbation theory, a variant of chiral per¬ 
turbation theory that builds in the consequences of the 
contracted SU(2AT/) (including the degeneracy of the nu¬ 
cleon and the A)|55H5?]. Given this, it is highly plausi¬ 
ble that all examples of relations in Skyrme type models 
that hold independently of the parameters of the model 
are, in fact, true model-independent results that hold for 
group-theoretical reasons at large N c . In what follows we 
shall assume that this remains true for nucleon-nucleon 
scattering and nucleon-antinucleon as well. 


III. NUCLEON-NUCLEON SCATTERING 

Nucleon-nucleon scattering at large N c was also first 
discussed in Witten’s classic 1979 paper[7]. Here again 
the analysis was first done in the context of the heavy 
quark limit for simplicity; it was argued that the natural 
language for describing such processes is time-dependent- 
mean-field theory (TDMFT). A key point in this analy¬ 
sis was that in order to have a smooth large N c limit in 
TDMFT, the initial conditions must be taken with fixed 
velocity as N c —► oo rather than with fixed momentum. 
This in turn means that the initial momentum is of order 
Nl , as is the initial kinetic energy. The regime of large N c 
with fixed initial velocity is referred to as “Witten kine¬ 
matics” . The analog to TDMFT for situations including 
light quarks in models such as the Skyrme model—which 
are expected to accurately encode the leading large N c 
scaling behavior and all consequences of emergent large 
N c symmetries—is time-dependent classical field theory. 

In ref. [Tj, Witten identified TDMFT as the appro¬ 
priate formulation of baryon-baryon scattering and also 
identified the appropriate kinematic regime being mo¬ 
mentum of order TV*. However, it is striking that ref. 
[7j did not identify precisely what is being calculated in a 
TDMFT calculation. It is noteworthy that TDMFT does 
not yield S-matrix elements—the basic quantum mechan¬ 
ical objects characterizing the scattering[9]. Similarly, 
one cannot directly compute the total cross-section, the 
total elastic cross-section or the differential elastic cross- 
section. Rather, as noted in ref. [5], TDMFT (or the 
analogous time-dependent classical solutions in Skyrme 
type models) allow one to compute variables associated 
with flows of conserved quantities such as baryon num¬ 
ber or energy[5]. (Although, it was subsequently shown 
that by studying the breakdown of the regime of validity 
of TDMFT as a function of impact parameter, one can 
get information about the total cross-section and elastic 
cross-section[10, IT]) Using TDMFT, in Witten kinemat¬ 
ics one can, for example, meaningfully compute the an¬ 
gular dependence of the outgoing energy relative to the 
beam axis, by doing TDMFT calculations with different 
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impact parameters and averaging over impact parame¬ 
ters. In essence, these quantities correspond to averages 
over many S-matrix elements. 

At first sight, this seems quite promising: there are 
physically relevant quantities which can be obtained from 
mean-field or classical dynamics in Witten kinematics. 
Unfortunately, in practice this is of limited utility. Note 
that the explicit form of the TDMFT equations derived 
in ref. [2 are only valid in the regime in which all quark 
masses are much larger than Aqcd- For realistic quark 
masses of relevance to physical nucleons, the form of the 
TDMFT for QCD equations is unknown. Of course, one 
can always use a model such as the Skyrme model in place 
of QCD. However, such models are not QCD and thus 
cannot reproduce the detailed results of large N c QCD. 
Such models have an important virtue however. They are 
believed to faithfully reproduce the leading N c scaling of 
observables and correctly encode both any approximate 
symmetries of QCD (such as chiral symmetry) and the 
leading order effects of emergent spin-flavor symmetry of 
baryons. Thus, one can use models of this sort as a tool 
to identify results of large N c QCD which depend only 
on the symmetry properties and scaling laws. We note in 
passing that for certain baryon observables the large N c 
and chiral limits are not uniform and results may depend 
on the order of limits [5HH2II ; the semi-classical analysis 
of soliton models corresponds to taking the large N c limit 
first with fixed quark mass prior to the chiral limit; one 
can then do a chiral expansion for each coefficient in a 
1/N C expansion. 

Topological solitons in chiral field theories based 
on hadronic degrees of freedom such as the Skyrme 
model [2TH23l can be used to identify model-independent 
properties of scattering as well as static properties. The 
analysis to do this was developed in [5] and this approach 
will be briefly reviewed here. It turns out that the pri¬ 
mary predictions involve the spin-flavor dependence of 
various observables. Of course, it is precisely these that 
we expect to be associated with the contracted SU(21V/) 
symmetry. Let us see how to use these models to deduce 
the spin-flavor dependence of variables associated with 
the collective flow of conserved quantities such as energy, 
momentum, or baryon number. One begins the analysis 
by imagining that one has access to an arbitrarily large 
set of solutions of the classical (or mean-field) scatter¬ 
ing processes associated with a given initial velocity that 
fixes the incident energy. These classical solutions all in¬ 
volve initial conditions of two widely separated rotated 
hedgehog baryons moving towards each other offset by 
an impact parameter b, a two dimensional vector. For 
concreteness since the experiments we have in mind are 
typically fixed-target experiments, we can restrict our at¬ 
tention to situations in which one soliton is stationary 
and the other is moving toward it. Thus, the initial con¬ 
ditions depend on the velocity, the scattering axis and 
8 more collective variables—2 for the components of the 
impact parameter vector and 6 to specify the orientation 
of the two initial hedgehogs. 


The classical solutions are, in effect, movies. The key 
issue is how to extract from these movies physically rele¬ 
vant observables in the underlying quantum theory which 
hold to leading order in 1 /N c . The first step in doing so 
is to find quantities which can be identified clearly in the 
classical calculation (movie) which has physical mean¬ 
ing even in the quantum context. An example might be 
the fraction of the initial kinetic energy of the baryons 
which gets converted into energy in mesons (including 
the meson masses). In any given classical solution of 
the Skyrme model this could be identified and similarly 
in any physical scattering process this could—at least in 
principal— be measured. The natural way to quantify 
the physical observable in scattering is in terms of some 
sort of cross-section. Thus for example one could ask 
for the cross-section for collisions in which the fraction 
of initial energy is converted into mesons is greater than 
some fixed value, /. This quantity is well-defined ex¬ 
perimentally. Quantum mechanically it can be obtained 
by summing over a large number of S-matrix elements 
(the number will scale with N c since kinematically the 
number of open channels that can contribute scales with 

N c ). 

The central point is that this quantity is also calculable 
classically and at large N c \ the classical result correctly 
accounts for the quantum results up to 1 /N c corrections 
(assuming that the model correctly describes the QCD 
physics). The cross-section is calculable classically in the 
following way: for any given classical solution with given 
impact parameter one can identify the fraction of the ini¬ 
tial kinetic energy converted to meson field energy. One 
can map out the region in impact parameter space for 
which this fraction is greater than /. The classical cross- 
section is simply the area in impact dimension space 
corresponding to such solutions. Note that this cross- 
section is independent of N c in the large N c limit since 
the “movie” associated with the solution is independent 
of N c . Note moreover that an analogous procedure can 
be done to compute the cross-section for any property 
which is similarly independent of N c . 

However, this procedure does not yet give the cross- 
section of physical interest. Note that the classical calcu¬ 
lation still depends on the variables associated with ini¬ 
tial orientations of the hedgehogs, A\ and A 2 . Thus, the 
cross-section associated with some property V computed 
classically is a function of A\ and A 2 : a v (Ai, A 2 ;p) 
where p is the initial momentum and of order N*. In 
contrast the physical scattering observables are typically 
given in terms of the spin and flavor of the incident 
baryons with some spin quantization axis: , 2 - 

J 1 1 rrig 

Since the directions of the spin projections are in princi¬ 
ple arbitrary, the natural way to write the cross-section 
is as the expectation value of a quantum mechanical op¬ 
erator in spin-isospin space. The most general form con¬ 
sistent with isospin symmetry, time-reversal and parity 
is 
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o-^i m i m 2 m? (p) = {mlml m 2 s m 2 i \d v (p)\mlm\m 2 s m 2 ) with (3) 

a v (p) = X 0 (p) + Y 0 (p)ai ■ cr 2 + Z 0 (p)(a i • h)(a 2 • n) + (Xi(p) + Ti(p)cri • a 2 + Z 1 (p)(a 1 ■ n)(a 2 ■ n))r \ • r 2 


where the hat in dr v (p) indicates that it is a quantum 
operator in spin-flavor space; in contrast, the hat in n 
indicates that it is a unit vector which is in the direc¬ 
tion of of the scattering axis, p is the momentum of 
the particle in the lab frame. Thus, there are in gen¬ 
eral six functions of the initial velocity (and possibly 
other variables associated with the observable) that fully 
characterize the cross-section associated with property 
V: Xq, X\. yjj, Y \, Z() and Z\. 

The critical issue is to relate cr v (Ai, A 2 )(p), the 
classical cross-section parametrized by the orientation 
variables A, with X 0 (p), X 1 (p),Y 0 (p),Y 1 , (p), Z 0 (p) and 
Zi(p). To do so, it is sufficient to recall that the initial 
orientation variables A are collective coordinates which 
are adiabatic and at large N c decouple from the other 
degrees of freedom in the problem. Hence, once can pro¬ 


mote these to quantum variables and requantize the mo¬ 
tion associated with them. The act of doing so converts 
classical configurations which break both spin and fla¬ 
vor into states with well defined spin and flavor quantum 
numbers. This means that if the initial nucleons are in 
some spin-flavor state | ip), then up to 1/N C corrections, 
the cross-section associated with V will be given by 

°%(p) = (V'l'^IV’) (4) 

= j dA 1 d,A 2 |V>(Ai, A 2 )\ 2 a' p (A 1 ,A 2 ;v) 

where the function cr p {Ai, A 2 ;p) is obtainable classi¬ 
cally. The collective wave-functions for states associated 
with well-defined spin and isospin are simply normalized 
Wigner D-matrices. This implies that the leading order 
expression is given by 


mi mf mi m. 


2 = (m, ttia m q tyia \a uia ra„ m: 


l ) = 


JdA ld A 2 \Dl_ AA^\dL 


(5) 


The fact that the collective wave-functions are sim¬ 
ply normalized Wigner D-matrices allows one to deduce 
a critical result: at leading order in 1/N C , must be 
invariant under a simultaneous flip of spin and isospin 
for either of the incident nucleons. This follows from a 
well-known property of Wigner D-matrices: 

= (—!)“' (Di mi _ m ,(H))* (6) 

which implies that \D J mrn ,{A)\ 2 = \D J _ m _ m ,(A)\ 2 . Com¬ 
bining this with Eq. implies that at leading order in 
the 1/N C expansion 

^m\m^ mlm 2 ^P) ^—m\—m 1 ml m 2 ^P) 

‘ P \ , (7) 

ml —ml — ml 'P J 

i.e. a simultaneous flip of spin and isospin for either of the 


incident nucleons leaves the cross-section invariant. On 
the other hand, in general \D 3 mTn ,(A)\ 2 7 ^ |Di mm ,(H )| 2 
and |D^ m ,(H )| 2 ^ \D^_ m ,(A)\ 2 . Thus 

a ml m] ml ml ip) ¥= a -ml mj m 2 s ml (?) 
a ml ml ml ml (p) ¥= -ml ml ml (?) 

a ml ml ml ml ip) ¥= m J -m 2 s ml 0 3 ) 

a ml ml ml m ? (p) ¥= mj m 2 s -ml (?) 

i.e. a flip of only the spin or only the isospin for one of 
the two nucleons alters the cross-section at leading order. 

Comparing Eqs. ([7]) and Q with Eq. ([ 3 ]), it is apparent 
that of the six functions which parameterize the cross- 
section in general, Xq. X\, Yq,Y-[ . Zq and Z\. only three, 
X 0: Y[ and Z- t contribute at leading order in the 1 /N c 
approximation: 


a v = 


(x 0 (p) + (bi(p)cn • cr 2 + Z 1 (p)(ai ■ h)(a 2 ■ n))n • r 2 ) x ^1 + 0 


(9) 


Equation has significant predictive power at least if the 1/N C corrections are negligibly small. Note that 
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half of the possible terms allowed by rotation and isospin 
symmetries are absent in the leading order expression. 
This means that certain quantities which could in prin¬ 
ciple differ will be identical at large N c . 

As an example, consider scattering where either the 
beam or the target (or both) are unpolarized. Even if 
only the beam or the target is unpolarized it is clear 
from the general structure of scattering in Eq. (|3| that 
polarization of the other particle is irrelevant. Since po¬ 


®pp unpolarized \P) = ®np unj 


larization is irrelevant, only Xq and X\ can contribute to 
scattering in Eq. ([T]). Of course this structure encodes 
isospin invariance, which implies that a^ p unpolarized (p) = 

a nn unpolarized (?)• However if X x ^ 0, as is generally ex¬ 
pected in the absence of a symmetry reason for it to van¬ 
ish, one sees that aJ punpola ^ ized (p) ^ ^ unpolarized (p)- 
However, as the large N c limit is approached, a con¬ 
tracted spin-flavor symmetry emerges and X\ — > 0. Thus 
for large N c , Eq. ^ implies that 


(p)x ( 1 + o (£0 (10) 


which is a prediction. Once polarized beams and targets 
are considered a large number of similar predictions can 
be madeP - It should be clear that the predictive power 
is due to the emergent spin-flavor symmetry. Note the 
analysis holds regardless of the parameters in the Skyrme 
Lagrangian and thus there are strong reasons to believe 
that it is a true model-independent prediction. 

It is worth recalling that predictions such as the one 
in Eq. ( fIo| ) are valid only in Witten kinematics which 
means that the incident momenta are taken to be of or¬ 
der N c and for cross-sections associated with properties V 
which can be computed in classical or mean-field treat¬ 
ments. Recall also that V is typically associated with 
bulk properties involving energy or the flow of baryon 
number. Thus for example, one can let V indicate the 
cross-section for scattering in which the baryons are de¬ 
flected by more than some fixed angle (with any number 
of pions produced). By differentiating with respect to 
the angle one can obtain differential cross-sections. How¬ 
ever, it should be stressed that these differential cross- 
sections are not the usual differential cross-section for 
elastic scattering but rather a kind of semi-inclusive dif¬ 


ferential cross-section which has a fixed angle for the 
outgoing nucleons but includes any number of created 
mesons. 


IV. NUCLEON-ANTINUCLEON SCATTERING 

The problem of nucleon-antinucleon scattering was 
also considered in Witten’s classic paper. He again jus¬ 
tified a time-dependent mean-field approach. The cen¬ 
tral observation at the core of the present paper is that 
the spin-flavor dependence of certain classes of nucleon- 
antinucleon scattering observables in Witten kinemat¬ 
ics can be obtained in an analysis nearly identical to 
that of nucleon-nucleon scattering. The analysis of 
nucleon-antinucleon scattering applies to qualitatively 
the same types of observables, cross-sections associated 
with bulk properties involving energy or the flow of 
baryon number— i.e. things computable at the classi¬ 
cal or mean-field level—in the Witten limit. It also leads 
to the same type of spin-flavor dependence: 


V ( \ / N N „,N „,N|-'P/ M N 

CT mN m N m N m NW = K rrii ™ s m i ]cr (p)\m s m l 
a v {jp) = X 0 (p) + Y 0 (p) ct n • CTjf + Z 0 (p)(a N ■ h)(a w 
where X 0 (p) ~ V c 0 , Y^p) ~ iV c °, Z x [p) ~ N° c 


m^mf) with 

• n) + (Xi (p) + Y 1 {p)a N ■ o-pf + ^i(p)(cr N • n)(a^- n)J 

md z <m~w c 


in) 


Tn * 'T’n 


The logic underlying this is analogous to the nucleon- 
nucleon case. Once again, for the purpose of identifying 
model-independent relations, models based on topolog¬ 
ical solitons serve as a surrogate for the full problem. 
Again one begins the analysis by imagining access to an 
arbitrarily large set of solutions of the classical (or mean- 


field) scattering processes associated with a given initial 
velocity. In this case, these classical solutions all involve 
initial conditions of two widely separated rotated hedge¬ 
hog baryons, one with winding number one (a baryon) 
and the other with winding number negative one (an an- 
tibaryon). They move towards each other offset by an 
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Relation number 

Annihilation Cross-Section Relation 

1) 

^pn unpolarized (.P) ®pp unpolarized (.P) 

HEM 

2) 

<7mn(P) = x 

'l + O 

£ 



3) 

ffp^W = CT PPtt(P) x 

'l + O 

* 



4) 

2 (^pn 'f'l'Cp) ^pp' |-'|'(p)) ®pn unpolarized (.P) ^ ( 


M 

5) 

| (^pnV (?) ~ a pn L u(P)) ~ “I ( a pp L n(P) ~ a pini(P)) 

X (l + O ( 

m)) 

6) 

Vpn T n(p) = a pp T n(p) x 

'l + O 

re 

) 


7) 

ffpn T u(p) = a pp T n(p) x 

'l+O 

re 

; 


8) 

2 (^prifT^P) “1“ ^ppTT^)) ^pn unpolarized (.P) f 


M 

9) 

2 (^pn ft Op) ^pn t4- (p )) 2 (P’pptt^P) a PP T4-0p)) 

x (l + 0\ 

M 


TABLE I. Relations between the total annihilation cross-section for nucleon-antinucleon scattering in various spin and flavor 
channels that hold for sufficiently large N c . The superscript A in a A indicates that the cross-section is for the total annihilation 
cross-section. The superscripts L and T stand for longitudinal and transverse respectively and indicate the polarization axis; 
longitudinal corresponds to spins polarized along the beam axis and transverse is for spins polarized perpendicular to it. Thus, 
°psTT corresponds to a proton and antineutron with both spins + | quantized along the positive beam direction. The subscript 
“unpolarized” corresponds to situations when either the nucleon or the antinucleon (or both) are unpolarized. 


impact parameter b, a two dimensional vector. It is nat¬ 
ural to restrict one’s attention to situations in which the 
baryon is initially at rest and the antibaryon moving to¬ 
wards it, as this is the typical experimental set up. Thus, 
as in the case of nucleon-nucleon scattering, the initial 
conditions depend on the momentum, the scattering axis 
and 8 more collective variables—2 for the components of 
the impact parameter vector and 6 to specify the orienta¬ 
tion of the two hedgehogs. Again for fixed values of the 
orientation variables, one can map out which values of 
the impact lead to classical solutions satisfying property 
■p; the area in impact parameter space defines the cross- 
section cr v (Ai, A 2 ;p). Again the orientation degrees of 
freedom are requantized. Finally, a comparison of the 
most general form for the cross-section to the form ob¬ 
tained from taking matrix elements of a^ > (A ll A 2 ;p) us¬ 
ing the requantized A variables and exploiting Eq. ([b]) 
leads to Eq. (If) 


V. TOTAL ANNIHILATION CROSS-SECTION 

As noted in the introduction, the type of observables in 
nucleon-nucleon scattering for which the spin-flavor de¬ 
pendence is determined at large N c are not those which 
experimentalists analyze and report. The analogous ob¬ 
servables in nucleon-antinucleon scattering are also typ¬ 
ically not analyzed in experiments. However, nucleon- 
antinucleon scattering has one observable that has no 
analog in nucleon-nucleon scattering, is in the class of 
observables for which the large N c analysis applies and 
is commonly analyzed in experiment: the total annihila¬ 
tion cross-section. 

It is easy to see that the total annihilation cross-section 
is in the class of observables for which the analysis ap¬ 
plies. In effect, this means that it can be calculated clas¬ 


sically for any given initial classical configuration param¬ 
eterized by the hedgehog, orientation angles, the impact 
parameter vector and the velocity. Starting with an ini¬ 
tial configuration in the Skyrme model, one can follow 
the held configuration forward in time. The baryon den¬ 
sity ps is given by the topological winding number den¬ 
sity, which is fixed by the field configuration. Initially 
the baryon density distribution consists of two well sepa¬ 
rated blobs, one with a baryon number one and the other 
with baryon number of negative one, which are traveling 
towards each other, off-set by the impact parameter. As 
time goes forward, the distribution of baryon density may 
get complicated as the baryon and antibaryon interact. 
However, at sufficiently long times after the interaction 
begins, things will necessarily simplify: as t —» oo, either 
Pb , the baryon density, will approach zero everywhere 
(corresponding to annihilation) or it will be in the form 
of two lumps, one with integrated baryon number one and 
with baryon number of negative one moving away from 
each other with a velocity whose magnitude is less than 
the initial velocity (corresponding to an inelastic scatter¬ 
ing process in which energy is lost to meson emission). 
Which of these two outcomes occurs is completely deter¬ 
mined by the orientation angles, the impact parameter 
vector and the momentum—this is a necessary conse¬ 
quence of the calculation being classical and hence deter¬ 
ministic. The classical annihilation cross-section for fixed 
p, A-[. and A 2 is simply the area in impact parameter 
space of initial configurations for which the ps (x 7 1) — > 0 
for all x as t —> oo. 

The upshot of this is that the spin-isospin dependence 
of the total annihilation cross-section at large N c is given 
by Eq. ©• To the extent that the subleading terms 
in 1/N C are negligible, it means that all total annihila¬ 
tion cross-sections are expressible in terms of just three 
functions of the initial velocity. This enables one to re- 















late the total annihilation cross-section in many different 
spin-isospin channels up to corrections of relative order 
1/N C . A number of these relations are given in Table 
[TJ All of the relations in Table |T] are true predictions of 
large N c QCD; that is, these relations depend on more 
than mere isospin and rotational invariance. A quick 
glance at the table makes clear that relations 2 and 6, 3 
and 7, 4 and 8, and 5 and 9 have identical forms. The 
only difference is that one relation applies for longitu¬ 
dinal polarization and the other applies for transverse. 
However, despite the formal structures, they are distinct 
predictions-in general, the annihilation cross-section for 
transversely polarized nucleons and antinucleons differ, 
even at large N c . 

The list of relations in Table [I] is not complete. Clearly 
one can exploit rotational invariance or isospin invariance 
to generate new relations from those in the table. Thus, 
if one simultaneously switches all protons with neutrons 
and antiprotons with antineutrons, all the relations re¬ 
main true; similarly if one swaps all spin up states and 
spin down for both particles simultaneously, all of the re¬ 
lations remain true. One can also generate new relations 
by taking linear combinations of the existing relations. 
Finally one can also generate new relations from Eq. 
for situations where the spin polarization directions are 
neither longitudinal nor transverse but at any given an¬ 
gle to the beam and for cases where the polarizations of 
the nucleon and antinucleon are quantized along different 
axes. 

One of the principle reasons for the focus on the to¬ 
tal annihilation cross-section—apart from the fact that 
it is an intrinsically interesting quantity—is that it is one 
which can be studied experimentally. Unfortunately, to 
the best of our knowledge, experiments measuring anni¬ 
hilation cross-section with polarized antineutron beams 
have not been done. This is hardly surprising since such 
an experiment would be immensely challenging. This 
means that most of the relations in Table U cannot be 
tested empirically. Fortunately, relation 1) does not de¬ 
pend on polarized beams and ought to be testable. It pre¬ 
dicts that the total annihilation cross-section for proton- 
antiproton and proton-antineutron reactions should be 
the same up to l/iV c corrections when the system is in 
Witten kinematics. 

Before confronting relation 1) with experimental data, 
it is important to clarify the nature of Witten kinematics, 
i.e. the requirement that p is of order N c . Formally 
this means that a nucleon-nucleon or nucleon-antinucleon 
cross-section associated with property V in spin-isospin 
channel c satisfies 

*f(p) = y?(fl x (i + o(A)) (l2) 

with p = p/N c 

where (p) is independent of N c . It is important to 
note that the coefficient in front of the relative order 
(1 correction can depend both on the V, the property 


defining the cross-section, and p = p/N c . It is possible 
that this coefficient grows with decreasing p, and indeed, 
that for certain choices of V, the coefficient diverges as 
p —> 0 indicating a breakdown of the 1/N C expansion at 
p = 0. This can happen for the following reason: at any 
N c , including arbitrarily large ones, p = 0 corresponds 
to p = 0 while for any non-zero value of p at sufficiently 
large N c , p is much larger than characteristic hadronic 
scales which such mesons masses and the size of the nu¬ 
cleon which might be expected to control the onset of 
semi-classical/mean-field regime on which the large N c 
analysis is based. 

Given this, there is a concern that at sufficiently low 
p for any given value of N c , the system might be out¬ 
side the regime of validity of the 1/N C analysis and the 
spin-flavor predictions based on it cease to be valid. In 
the case of nucleon-nucleon scattering, it is clear that 
this happens essentially regardless of the property, V. 
Note that as p —» 0, the scattering becomes purely elas¬ 
tic s-wave scattering and involves a single quantum me¬ 
chanical channel. Since this is fundamentally quantum 
in nature, one expects that any semi-classical analysis 
should break down at low enough momenta. Empirically, 
it was seen in ref. [20] that the large N c predictions of 
the spin-flavor dependence failed when used in the elas¬ 
tic scattering of nucleon-nucleon scattering. If a similar 
breakdown of the validity of the large N c analysis occurs 
for nucleon-antinucleon scattering at low momentum, it 
will be difficult to assess empirically the predictive power 
of the large N c approach in nucleon-antinucleon annihi¬ 
lation, since the antineutron-proton annihilation data is 
at rather low beam momenta—below 500 MeV. 

Fortunately, low momentum nucleon-antinucleon an¬ 
nihilation differs from low momentum nucleon-nucleon 
scattering in a fundamental way. While nucleon- 
nucleon scattering at very low momentum involves only 
a single quantum channel (the elastic s-wave channel), 
nucleon-antinucleon annihilation always involves many 
channels—even as the incident momentum goes to zero. 
The decay can go into two mesons in s-wave, three mesons 
in various angular momentum and isospin combinations, 
4 mesons in multiple combinations, etc. As the N c limit 
is approached, the number of such channels grows rapidly 
as a function of N c . Since the crux of how a classical re¬ 
sult emerges from a quantum scattering problem is the 
contributions from many quantum channels which be¬ 
have similarly, and the classical regime was key to the 
derivation of the spin-flavor relations for scattering, it is 
plausible that what is driving the validity of the spin- 
flavor dependence for scattering is the presence of many 
kinematically available quantum channels. If this is the 
case, then one expects the relations in Table [I] to hold at 
large N c for arbitrarily small values of p. While this is not 
guaranteed to be correct, it is plausible and motivates a 
comparisons of relation 1) with the available data, even 
though the data is at low momentum. 

In Fig. [l] data from ref. [28] for annihilation cross- 
section of an unpolarized antineutron beam on a proton 
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FIG. 1. The data points represent the total annihilation cross- 
section for unpolarized antineutrons incident on protons with 
various beam momentum; the data comes from ref. [28] and 
the errors include both systematic and statistical errors added 
in quadrature. The solid curve is a phenomenological fit to 
the total annihilation cross-section for unpolarized antipro¬ 
tons incident on protons as a function of beam momentum 
from ref. [29]. This curve accurately fits all of the data points 
in this range to within a few percent. 

target is given as a function of beam momentum. To the 
extent that relation 1) of Table [I] applies, the antiproton- 
proton annihilation cross-sections should be the same up 
to creations of relative order 1 /N c . For ease of compar¬ 
ison, rather than giving individual data points for the 
antiproton-proton annihilation cross-section, a fit to this 
data taken from ref. [25] is given. The curve is entirely 
phenomenological in nature and is not based on any un¬ 
derlying theory. However, the quality of the fit is ex¬ 
tremely good. The key point is that the quality of this 
fit is sufficiently good that for all data points, the dis¬ 
agreement between the fit and the data is much smaller 
than the quoted error bars for any of the antineutron- 


proton data points. It thus seems that the fit is a useful 
basis for comparison with the antineuton-proton data. In 
practice, the antiproton-proton data points at the lower 
end of this range are accurately fit to within a few per¬ 
cent by the curve and most of the points are fit to better 
than a percent. 

The comparison between the fitted antiproton-proton 
data and the antineutron-proton data is rather striking. 
By eye, it seems that the expectation based relation 1) 
holds. Indeed, the fit looks as though it could have been 
based on the antineutron-proton data rather than the 
antiproton-proton data. One way to quantify this: the 
X 2 per degree of freedom of the anti-neutron data taking 
the curve fit to antiproton-proton as the theoretical pre¬ 
diction is less than unity. This is remarkable, in that even 
if relation 1) is valid, it only expected to hold to order 
1/N C which for the real world is 1/3. Thus, given given 
the quality of the data, it seems quite safe to conclude 
that empirically relation 1) does hold. 

Of course, the fact that the unpolarized antiproton- 
proton and antineutron-proton total annihilation cross- 
sections are so similar to what is expected at large N c 
does not necessarily indicate that this is due to the viabil¬ 
ity of large N c physics for scattering at N c = 3. It is quite 
possible that the agreement with the large N c prediction 
is accidental and stems from some other cause. The case 
would be strengthened significantly if the data extended 
to significantly higher momenta (say up 2 GeV) where 
the question of whether the momentum is high enough 
for system to be in the regime of Witten kinematics does 
not arise. More compelling would be data with polarized 
beams acting on polarized targets using both antineutron 
and antiproton beams. This would allow tests for more 
than just relation 1). However, despite the limitations of 
the data, it is nevertheless encouraging that the data we 
do have is consistent with the large N c expectations. 
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